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ABSTRACT 

We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly colli¬ 
mated relativistic jets can be produced in neutron star merger models of short-duration gamma ray 
bursts without the need for a highly directed engine or a large net magnetic flux. Even a hydro- 
dynamic engine generating a very wide sustained outflow on small scales can in principle produce 
a highly collimated relativistic jet, facilitated by a dense surrounding medium which provides a co¬ 
coon surrounding the jet core. An oblate geometry to the surrounding gas significantly enhances 
the collimation process. Previous numerical simulations have shown that the merger of two neutron 
stars produces an oblate, expanding cloud of dynamical ejecta. We show that this gas can efficiently 
collimate the central engine power much as the surrounding star does in long-duration GRB models. 
For typical short-duration GRB central engine parameters, we find jets with opening angles of order 
10 degrees in which a large fraction of the total outflow power of the central engine resides in highly 
relativistic material. These results predict large differences in the opening angles of outflows from 
binary neutron star mergers versus neutron star-black hole mergers. 

Subject headings: hydrodynamics — relativistic processes — shock waves — gamma-ray burst: general 
— ISM: jets and outflows 


1. INTRODUCTION 

Gamma ray bursts (GRBs) i nvolve energ e tic, colli¬ 
mate d relativistic outflows (see iPiranI 12004 iMeszaro^ 
I2OOOI for recent reviews). The current consensus on 
short-duration GRBs is that they are produced during 
the merger of a binary neutron star system or a neutron 
star-black hole system. However, the detailed proper¬ 
ties of the central engine which powers the burst remain 
uncertain. 

GRB engine models often require both fast rotation 
and large magnetic fields, in order to produce a highly 
collimated jet. However, the mechanism for producing 
a sufficiently strong large-scale field is uncertain, lead¬ 
ing to the possibility that it might not be present. For 
example, the Blandford-Znajek mechanism typically re¬ 
quires 10^® Gauss in large-scale poloidal magnetic field 
to power a GRB, but such large magnetic flux is not ex¬ 
pected to be available in typical neutron star mergers. It 
may instead be the case that the magnetic field in the 
central engine remains toroidal and/or small-scale, and 
is primarily responsible for providing angular momen¬ 
tum transport in the accretion disk. From an observa¬ 
tional perspective, beaming- corrected GRB energie s are 
inferred to be - lO^s-so erg (iFong et a1 .1127111 Infill . yet 
calculations of typical post-merger disk masse s suggest 
-MdiskC^ ~ 10®^ erg l|Hotokezaka et all l2013af) . imply¬ 
ing a very inefficient engine without a strong large-scale 
magnetic fl ux threading it (contrar y to what is often as¬ 
sumed, e.g. iPaschalidis et al.ll2ni4 l. In fact, in the long 
GRB context, even if one assumes the existence of strong 
large-scale poloidal magnetic fields, it can still be difficult 
for some engine models to produce a highly collimated 
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outflow (|Mdsta et al.ll2ni4 l. 

These concerns are severe enough that we take seri¬ 
ously the possibility that strong, large-scale magnetic 
fields may not be present in a short-duration GRB cen¬ 
tral engine. In this study, we investigate an alterna¬ 
tive picture, in which the engine itself does not need to 
be strongly collimated, but the outflow becomes a col¬ 
limated jet after interaction with a dense surrounding 
medium. This idea has been previously proposed to ex¬ 
plain t he formation of radio com ponents in active galactic 
nuclei (|Blandford &: Re^ 11974) and is known to be an 
important part of the collima tion process in long GRBs 
(|MacFadven fc Wooslevl[T999ll . 

The process of collimation from a wide engine to 
a narrow jet can be enhanced if the surrounding 
medium has an oblate geometry (as also suggested by 
iBlandford fc Ree^ 119741 . Oblate density profiles are 
quite natural in GRB scenarios, as most engine models 
require rapid rotation for operation. 

Under the right circumstances, even a spherically sym¬ 
metric engine can produce a highly collimated outflow if 
the geometry of the surrounding medium is sufficiently 
oblate. This is ultimately due the fact that a blastwave 
expanding in an asymmetric external medium will pref- 
erentially favor the direction of steepest density gradient 
(lKomnaneetslll96Til : iKorvcanskvl 1 1 9921 : iHnatvk fc PetrukI 
I1999II . Thus, an oblate external medium can result in 
a prolate blastwave. If the engine is persistent, contin¬ 
uing to push from the inside of this prolate bubble, the 
aspect ratio of the bubble will become increasingly exag¬ 
gerated until the flow is collimated to within a sufficiently 
small opening angle. At this point, recollimation shocks 
can take over (iMacFadven et al.l[200ll : IZhang et al.ll2003t 

iBromberg et al.ll2011ll and maintain the collimation until 

the jet completely breaks out of the surrounding cloud. 

In the present study, this collimation process will be 
demonstrated numerically, in the context of short GRBs. 
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Numerical relativity simulations show that shortly after a 
neutron st ar merger, an oblate cloud of dynamical ejecta 
is present (iHotokezaka et aLll2013bf) . The idea that this 
cloud can provide a collimating cocoon has also recently 
been investigated (iNagakura et al.l[2ni4[l . though in that 
study the parameter space of engines was not thoroughly 
explored and the surrounding cloud was assumed to be 
spherical. Starting with this dynamical ejecta cloud as 
initial conditions, we investigate the collimation process 
for a range of engine parameters and including the oblate 
geometry, and find that very wide hydrodynamic central 
engines are capable of producing collimated jets. 

2. NUMERICAL SET-UP 

This study performs a numerical integration of the 
equations of two-dimensional (2D) axisymmetric rela¬ 
tivistic hydrodynamics: 

df,{pu>^) = Sd (1) 

d^iphuf^u’'+ = S’' ( 2 ) 

where p is proper density, ph = p + e + P is enthalpy 
density, P is pressure, e is the internal energy density, 
is the four-velocity, and g^’' is the Minkowski metric 
(the equations are expressed in units such that c = 1). 
The source terms Sd and S’' will be used to model the 
injection of mass, momentum, and energy by the central 
engine on small scales. The equation of state is assumed 
to be radiation dominated: e = 3P. 

The equations are integr ated using the moving-mesh 
hydro dynamics code, JET (|Duffell fc MacFadvenll20lTL 
[20il). The JET code’s moving mesh makes it effec¬ 
tively Lagrangian, which allows for accurate resolution 
of highly supersonic and relativistic flows. The inner and 
outer boundaries are also moved during the flow’s evo¬ 
lution, so that the dynamic range covered at any given 
time is kept modest. The spatial resolution is given by 
A0 ^ Ar/r ~ .006. This value is not exact, as zones can 
expand and contract, but zones are actively refined or 
de-refined so as to maintain aspect ratios of order unity. 
Additionally, zones are more concentrated near the sym¬ 
metry axis, so that resolution there is about twice as high 
(A6» ~ Ar/r ~ .003). 

Initial conditions are modeled after the numer- 
ical output of a neut r on star merger sim ulation 
(|Hotokezaka et ahl l2013bl : INagakura et al.l 1201411 . It is 
an outflow described by the following: 
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where pc = 287Mo/R^, Ri = .06i?o, Ptrans = 

l0~^Mo/Rl, pisM = 10~‘^^Mo/Ro, and to = IORq/c. 
Unless explicitly stated otherwise, the characteristic 


scales are the initial radius of the cloud i?o = 850 km, and 
the mass of the cloud Mq = 10 “"^M 0 , though the hydro- 
dynamical equations are scale-invariant, and therefore 
the value of these two quantities is arbitrary, so long as 
all other dimensionful quantities are expressed in terms 
of these two variables (as well as the speed of light, c). 
In cgs units, the external density pisM ~ 3 x 10^ g/cm^ 
is much larger than typical ISM densities of ~ 1 proton 
per cm^, but this is irrelevant for the current study, as 
we are only concerned with the production of the jet and 
not its late-time deceleration. The density pisM is chosen 
to be sufficiently low so as to have negligible effect on the 
jet during the course of the calculation. The parameter 
w = 1.3 is the aspect ratio of the oblate cloud. 

A hydrodynamic engine is injected on small scales, us¬ 
ing t he source terms Sn and S’' in equations and ID), 
as in iDuffell fc MacFadvenI (|2014[1 . However, in contrast 
with that study, the engine in the present study will have 
a very wide half-opening angle, 6q = 0.5, or a full open¬ 
ing angle of about 60 degrees. This is comparable to the 
opening angle of the high speed outflow seen in hydro 
simulatio ns of radiatively ineffi cient black hole accretion 
disks le.g. lOhsuga et ahlfeoPSll or MHD simul ations with¬ 
out s ignificant large-scale magnetic flux (e.g. I Jiang et 
I2014f) . Additionally, in this study, the injection radius tq 
is a fixed multiple of the inner boundary, rg = 3r„iin- 

The engine duration is set to r = 30i?o/c = .085 sec¬ 
onds. The engine produces a clean relativistic outflow, 
with an energy-to-mass ratio of po = 10 ^, and a bulk 
Lorentz factor 70 = 10. The luminosity of the engine is 
given by 

Lq = eMoc^/r ( 8 ) 

where e is the ratio of the energy deposited to the cloud 
rest mass. Several calculations were performed, vary¬ 
ing the value of e from 10“^ to 0.24, implying a to¬ 
tal energy of Etot = eMgc^ in the vicinity of 10 ^'^ - 
IqSO ergs (again, the interpreted value of this energy 
depends on the choice of mass scale, Mq = 10“^Mq). 
The primary calculation assumes Etot = 4 x 10"^® erg, 
and Mq = 1 O“^M 0 , implying e = 0.024. Short GRB 
energies have bee n observationally inferr ed to be in the 
2Q48-50 erg range (iFong et al.ll201^ 1201411 . Alternatively, 
from th eoretical models, typical d isk masses are of order 
O.IM 0 (|Hotokezaka et ani2013af) . so that if this disk is 
powering the engine, the total energy available to the 
engine is ~ 10 ®^ erg, assuming ^ 10 % efficiency convert¬ 
ing the accretion power into an outflow. The mass of 
the surrounding cloud is expected to be in the range of 
10“^ — 1 O“^M 0 (jHotokezaka et ani2013bf) . and therefore 
the value of e = Etot/Moc^ could take on values any¬ 
where between 10 “^ and unity, comparable to the range 
explored in this study. 

3. RESULTS 

Figure [T] shows two example runs with Etot = 2.4 x 
lO'^Mgc^ = 4.3 X 10^® erg. The engine properties of 
these two runs are identical; the only difference between 
the two is the shape of the surrounding cloud, parame¬ 
terized by the aspect ratio, w. The left side of the figure 
shows a spherical setup with w = 1 and the right side 
uses w; = 1.3, a value chosen based on the aspect ratio of 
the dynamical ejecta found in binary neutron star merger 
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Fig. 1.— Comparison between oblate and spherical models for 
the surrounding cloud (spherical on the left side, oblate on the 
right). Dashed cyan curves indicate the expanding surface of the 
cloud of dynamical ejecta. Various moments in time (top to bot¬ 
tom, ct/Ro = 3, 5,15, 25) are shown as the engine attempts to clear 
a passageway in the cloud. The engine model is injected on a very 
small scale, which is only barely visible in this figure. The opening 
angle of the initial engine is very wide (full opening angle of about 
60 degrees), so the spherical cloud is not capable of collimating the 
flow by itself. The oblate cloud, however, is sufficient to produce a 
narrow jet. In both examples, the ratio of engine energy to ejecta 
cloud mass is Etot/Moc^ = 0.024. 

simulations (|Hotokezaka et aLll2013b[) . The w = 1 calcu¬ 
lation did not produce a collimated outflow, as the pres¬ 
sure from the cocoon by itself was not significant enough 
to keep the jet collimated. The final outflow is roughly 
spherical (see also Figure [2]). 

In the oblate {w = 1.3) run, the outflow from the en¬ 
gine is initially only weakly collimated, but this asym¬ 
metry becomes more prominent as the bubble inflates, 
preferentially moving along the steepest density gradi¬ 
ent. As the engine pushes, this asymmetry becomes in¬ 
creasingly exaggerated until recollimation shocks are ca¬ 
pable of keeping the jet collimated. The flow then clears 
a passageway, through which a highly collimated jet can 
escape. This collimated flow then breaks out of the cloud 
with a Lorentz factor in the tens, after which post-shock 
acceleration can increase this Lorentz factor into the hun¬ 
dreds. 
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Fig. 2.— Distribution of energy and terminal Lorentz factor with 
angle, for both spherical and oblate cloud models from Figure [T] 
The spherical model produces a nonrelativistic, spherical outflow, 
whereas the oblate cloud collimates the jet. For the oblate model, 
50% of the energy is contained within 8 < 650 = 0.12. The bottom 
panel shows an “averaged Lorentz factor” as a function of polar 
angle, defined as the ratio of energy to mass at the given angle. 
Etot/Moc^ = 0.024 in both examples. 

The remaining figures show how the flow is distributed 
in angle and Lorentz factor, and how the results depend 
on the energy injected by the central engine relative to 
the mass of the ambient medium (Figures [3] and [5]) and 
the injected opening angle of the outflow (FigureS]). The 
terminal Lorentz factor is measured via the specific en¬ 
thalpy, 7terminai = jh, which is the terminal Lorentz fac¬ 
tor if the fluid element is allowed to expand and adiabat- 
ically convert all of its thermal energy to kinetic energy. 
In the oblate example shown in Figure [T] about 30% of 
the energy has yh > 10. 

Figure [H shows the jet’s distribution with latitude. For 
the w = 1.3 example, 50% of the energy in the jet is 
contained within 0 = = 0.12. The w = 1 example is 

also included to show that the flow is roughly spherical 
and nonrelativistic if the surrounding medium is spher¬ 
ical (050 = 0.46 for the w = I model). Looking back 
to the fc = 1.3 case, the dip at very small angles im¬ 
plies a less relativistic jet core (this can also be seen in 
Figure [T]), but this is most likely an artifact of assuming 
axisymmetry. The maximum terminal Lorentz factor for 
this jet is around yhmax ~ 100. The “average Lorentz 
factor” ravg(0) (bottom panel of Figure [2]) is calculated 
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Fig. 3.— Dependence of terminal Lorentz factor, relativistic 
energy fractions, and opening angle of the jet on the total energy 
output of the engine, for a central engine opening angle of 6q = 0.5, 
and an oblate {w = 1.3) ejecta cloud mass of Mq = 10“^^ Mq. The 
relativistic hydro equations are scale-invariant, so these results can 
be generalized to different choices of Mq by holding Eioi/Mq fixed. 



Fig. 4.— Dependence of the jet half-opening angle, ^50, on the 
engine half-opening angle, 9q^ injected in the oblate cloud with 
w = 1.3 (e = Etot/MQC^ = 0.024). Engines with 6q < 0.5 (corre¬ 
sponding to a full opening angle of 60 degrees) produce collimated 
jets with 050 < 0.2, roughly independent of 6q. Wider engines 
are more difficult to collimate, though a surrounding cloud which 
is more oblate {w > 1.3) might facilitate collimation of wider en¬ 
gines. 



Fig. 5.— Comparison of the collimation by the spherical cloud 
(w = 1) with the oblate cloud {w = 1.3). The spherical cloud, at 
best, can collimate the jet to within 05o 0.25, but the oblate 

cloud is more effective at collimation. This model finds a minimal 
collimation angle of 05o 0.1 for the oblate cloud. 

via the ratio of total energy to total mass at a given polar 
angle. 

Figure [3] shows how these results depend on the total 
energy injected by the engine relative to the mass of the 
cloud, Etot/Moc^. For sufficiently low energy injection, 
Etot is negligible compared to the kinetic energy in the 
dynamical ejecta, and therefore the engine is unsuccess¬ 
ful at producing a collimated outflow. This leads to the 
low Lorentz factors and low energy in relativistic material 
for low values of i?tot in Figure [3] For sufficiently large 
Etot, the mass of the ejecta is not significant enough to 
collimate (or otherwise affect) the flow, and the outflow 
keeps the same opening angle it was given initially. Note 
that the limiting case in the high-energy limit is a jet 
with 7/1 ~ 10^ and ^50 ~ 0.5, roughly the same jet pa¬ 
rameters that were injected by the engine, implying that 
the surrounding cloud does little to affect the properties 
of the flow. 

The behavior in these two limits (low and high energy) 
imply an intermediate “sweet spot” for the total energy 
output of the engine, which is optimal for producing a 
collimated, relativistic outflow. For the model presented 
here, this “sweet spot” is e 10 “^, which corresponds to 
Etot in the few xlO"^® erg range for Mq = 1 O“^M 0 . This 
optimal energy likely depends on the detailed structure 
of gas in the surrounding cloud. For our engine and cloud 
model, at this sweet spot, about 10 — 30% of the energy is 
used to accelerate the relativistic material with yh > 10 , 
while the remaining 70 — 90% is deposited in the slower- 
moving material in the dynamical ejecta. 

Figure [4] shows the dependence of the final opening an¬ 
gle on 6 * 0 , the half-opening angle of the engine (for fixed 
Etot/Moc^ = 0.024). For sufficiently small 9o, the out¬ 
flow has 050 ~ 0.1 —0.2, roughly independent of Oq. How¬ 
ever, above some critical threshold, O^q begins increasing, 
so that the final outflow is wide when Oq > 0.7. 

Finally, Figure [S] compares the spherical {w = 1) and 
oblate {w = 1.3) ejecta models, to determine the extent 
to which oblateness affects collimation. The spherical 
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ejecta model has difficulty collimating such a wide en¬ 
gine, with a minimum ^ 0.25. The shape of the 
surrounding medium therefore plays an important role 
in the collimation process. 

4. DISCUSSION 

Dynamical processes during neutron star mergers cre¬ 
ate an expanding oblat e cloud of ejecta 
moving out at ^ 0.1c (iHotokezaka et ahl l2013bll . We 
have demonstrated that this ejecta collimates outflows 
from short-duration GRBs, even if the outflow has a large 
opening angle to begin with. There is thus no need for 
a central engine that itself generates a highly collimated 
relativistic outflow. This result is a variant of the well- 
known idea that hot flows can potentially self-coll imate 
as they become supersonic (jBlandford fc Ree^ll974[ l. We 
also show that the oblateness of the surrounding medium 
is an important ingredient to the dynamics. 

The canonical view of short-duration GRBs has been 
that they are powered by a relativistic jet a nalogous to 
radio - loud quasars or long-duration GRBs (jFong et al.1 
I2012L l2014fl . One indication that this analogy breaks 
down is that the outflow energies in short GRBs appear 
to be orders of magnitude below the energy available 
given remnant accretion disk masses of 0.1 Mr:^ found 
in sim ulations of neutron star mergers (|Hotokezaka et ^ 
[Mli). This is in contrast to other jet-dominated 
sources. This suggests that the black hole event hori¬ 
zon is not threaded by a dynamically important poloidal 
magnetic field ( i.e. the disk is not in a “magnetically 
arrested” state: iTchekhovskov et al.ll2011ll . Our results 
are consistent with this general inference in that colli¬ 
mated outflows from short GRBs need not originate from 
a central engine that in isolation produces a powerful col¬ 
limated outflow. 

Some small-scale and/or toroidal magnetic field is 
likely still an important part of the engine’s operation, 
but the field does not have to be responsible for colli¬ 
mating the flow to within Oj ~ 0.1, and the outflow does 
not have to be magnetically dominated. A full opening 
angle of about 60° appears sufficient for the surrounding 
medium to collimate the outflow from the central engine 
(Figures [Hand ID). 

This result removes some constraints from traditional 
short GRB engine models, and also opens up new pos¬ 
sibilities for what could be powering the central engine. 
Neutrino pair annihilation is more viable as an energy 
source if one does not require the engine itself to be colli¬ 
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mated. A 60° opening angle is also typical of disk winds, 
meaning that a disk wind could be a potential engine 
model, providing that a very clean (baryon-free) wind 
could be launched from the disk. 

The present study focused on the binary neutron star 
merger scenario for short GRBs. However, the same 
considerations apply to neutron star-black hole merg¬ 
ers, for which the dyn amical ejecta is exp ected to be 
much more oblate (e.g. iFoucart et al.ll20l^ . Based on 
the results presented here, we expect differences in the 
opening angles and Lorentz factors of outflows from black 
hole-neutron star mergers, relative to binary neutron star 
mergers. 

This result motivates several important future stud¬ 
ies. First, this same calculation should be carried out in 
full 3D, as it is possible that the collimation and ejec¬ 
tion process are sensitive to 3D effects. For example, the 
“hollow jet core” seen in Figure [2] is likely an artifact of 
the assumption of axisymmetry. Second, this calculation 
should be performed in the context of long GRBs, i.e. 
a Wolf-Rayet progenitor wh ich is oblate due to rotation 
([MacFadven fc Duffelll 1201^ in preparation). It seems 
likely that such a system will exhibit similar behavior. 
Third, larger parameter surveys are warranted to study 
specifically the oblateness of the surrounding medium, 
and its effect on collimation. It is also worth noting that 
Figure |T] shows a massive “plug” or “cork” in front of the 
jet, res ulting in a “top-heavy” flow, similar to what was 
seen in iDuffell fc MacFadvenl 1)20141 1. The properties of 
this cork and its impact on the jet dynamics should be 
investigated in more detail. 

Finally, it may be possibl e for analytical jet models 
(iBromberg et al.l l201l[ 1201411 to be modified to include 
the effect that oblate geometry can have on the prop¬ 
erties of a jet breaking out of a surrounding medium. 
This would be valuable for understanding what neces¬ 
sary properties an engine must have in order to produce 
a relativistic, collimated jet after breakout. 
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